Quantum Monte Carlo with Jastrow- Valence-Bond wave functions 
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We consider the use in quantum Monte Carlo calculations of two types of valence bond wave 
functions based on strictly localized active orbitals, namely valence bond self-consistent-field (VB- 
SCF) and breathing-orbital valence bond (BOVB) wave functions. Complemented by a Jastrow 
factor, these Jastrow- Valence-Bond wave functions are tested by computing the equilibrium well 
depths of the four diatomic molecules C2, N2, O2, and F2 in both variational Monte Carlo (VMC) 
and diffusion Monte Carlo (DMC). We show that it is possible to design compact wave functions 
based on chemical grounds that are capable of describing both static and dynamic electron correla- 
tion. These wave functions can be systematically improved by inclusion of valence bond structures 
corresponding to additional bonding patterns. 
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I. INTRODUCTION 

Quantum Monte Carlo (QMC) methods (see e.g. 
Rcfs. S-i) constitute an alternative to standard quan- 
tum chemistry approaches for accurate calculations of the 
electronic structure of atoms, molecules and solids. The 
most commonly-used approach consists of optimizing a 
flexible trial wave function in a variational Monte Carlo 
(VMC) calculation, and then using the resulting wave 
function in a more accurate fixed-node diffusion Monte 
Carlo (DMC) calculation. 

For atoms and molecules, the most common form of 
trial wave function consists of a Jastrow factor, expected 
to describe a major part of the dynamic electron correla- 
tion effects, multiplied by a single Slater determinant of 
orbitals expanded in a localized one-particle basis. For 
systems featuring important static correlation effects, the 
single Slater determinant is usually replaced by a lin- 
ear combination of several Slater determinants. This has 
been shown to lead to systematic improvements of both 
VMC and DMC total energies, provided that the wave 
function is properly reoptimized [^-|6[. In particular, full- 
valence complete active space (FVCAS) expansions, i.e. 
including all the Slater determinants that can be gener- 
ated by distributing all the valence electrons in all the 
valence orbitals, were found to yield near chemical ac- 
curacy for bonding energies of first-row homonuclear di- 
atomic molecules at the DMC level @. 

However, FVCAS expansions involve large numbers of 
Slater determinants that scale exponentially with system 
size, so this approach cannot be applied for large systems. 
Instead, more compact and yet systematic wave functions 
are desired for QMC calculations. Naively-truncated ex- 
pansions of Slater determinants of delocalized molecular 
orbitals cannot be seen as a general solution since they of- 
ten suffer from lack of consistency for calculating energy 



differences between two systems or two states. In con- 
trast, pairing-type wave functions should be an adequate, 
compact and systematic way of describing static correla- 
tion in QMC. Indeed, several forms of pairing wave func- 
tions have recently been used in QMC, namely antisym- 
metrized geminal power (AGP) B-jj, pfaffians [lOl. HH 
and generalized valence bond (GVB) |l2l |. with promis- 
ing results on atoms and small molecules. 

Here, we explore the use in QMC of valence bond (VB) 
wave functions based on nonorthogonal active orbitals 
that are localized on a single atom (for a review of VB 
theory, see Refs. IT3l - [l5[ ). The merits of valence bond 
self-consistent field (VBSCF) wave functions [l^, are 
tested, along with more flexible breathing-orbital valence 
bond (BOVB) wave functions [13; [3 • These approaches 
can describe static correlation effects with compact wave 
functions containing only a limited number of VB struc- 
tures (i.e., linear combinations of Slater determinants 
made of spin-singlet pairs of orbitals) which are selected 
on chemical grounds. In comparison to the GVB method, 
which is usually considered in the perfect-pairing and 
strong-orthogonality approximations, the VBSCF and 
BOVB approaches are more flexible and provide a greater 
chemical interpretability thanks to the use of localized 
nonorthogonal orbitals. In contrast to traditional ana- 
lytical calculations, nonorthogonal orbitals do not lead 
to extra computational cost in QMC. On the contrary, 
the use of localized nonorthogonal orbitals could poten- 
tially reduce the computational cost as compared with 
a delocalized molecular-orbital approach in three ways: 
a) localized orbitals are negligibly small and need not 
be computed beyond some cutoff distance, b) the num- 
ber of determinants needed may be reduced by having 
nonorthogonal orbitals, and, c) in the absence of spatial 
symmetry, the number of orbital parameters to be opti- 
mized is smaller for localized orbitals. 



BOVB wave functions improve upon VBSCF wave 
functions by allowing a different set of orbitals in each 

VB structure, effectively introducing dynamic correla- 
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than through backflow coordinate transformations [19[ or 
orbital-attached multi-Jastrow factors [l^l- In VB prac- 
tice, it is known that bonding energies at the BOVB level 
are significantly improved as compared with the VBSCF 
level, especially in cases where dynamic correlation is 
important. In QMC, this form of wave functions has 
already successfully been tested on acetylene dissocia- 
tion energies (2]| . However, it is not known whether the 
breathing-orbital relaxation yields significant improve- 
ment when the wave function is fully optimized in QMC, 
because it is conceivable that the Jastrow factor in QMC 
provides much of the same variational freedom. 

The paper is organized as follows. In Section II, we of- 
fer a brief review of the VBSCF and BOVB approaches, 
and describe our corresponding Jastrow- Valence-Bond 
wave functions used in QMC. In Section III, we present 
results for the four diatomic molecules C2, N2, O2, and 
F2, which span different types of bonding (single, mul- 
tiple, and three-electron bond) and electron correlation 
(strong dynamic correlation for O2 and F2, and strong 
static correlation for C2 and N2). We discuss the type 
of electron correlation to include in the wave functions 
to obtain accurate equilibrium well depths, as well as 
the importance of orbital optimization in QMC, and the 
computational cost of our wave functions. Section VI 
summarizes our conclusions. 



II. METHODOLOGY 

A. Valence-Bond wave functions 

The different valence bond methods can be classified 
into two main families, depending of the degree of local- 
ization of the active orbitals: semilocalized orbitals and 
strictly localized orbitals [22| . In the first family, the 
optimized active orbitals are expanded on all the basis 
functions of the entire molecule, and are thus free to 
fully dclocalizc, although most of the time they appear 
to be fairly localized. The GVB method [23l |24| and the 
spin-coupled valence bond (SCVB) method [23 belong 
to this family. In the second family, each active orbital is 
expanded in basis functions centered on a single atom or 
the atoms of a fragment of the molecule (alternative def- 
inition of strictly localized orbitals exists, Ref . [26l) . The 
VBSCF [H, and BOVB [13, \M methods belong to 
this family. The semilocalized approaches usually pro- 
vide more compact wave functions, but the delocaliza- 
tion tails of the orbitals may hinder their chemical read- 
ability. The strictly localized approaches usually involve 
more VB structures, but permit the distinction between 
covalent and ionic contributions to the bonds, thus pro- 
viding more insight into the nature of the chemical bonds. 
The interpretative power of the strictly localized VB ap- 
proaches is hence greater, and allows, among other ex- 
amples, the characterization of charge-shift bonding [28l - 
[soj . three-electron bonds [Sll - lssj or hypervalent bond- 
ing (ssl-lssj. Further, most applications involving VB 



curve-crossing diagrams [15|, |36[ are within the strictly 
localized VB formalism. 

Any VBSCF wave function can be written as a linear 
combination, 

|*vbscf) 

where |$/) are spin- adapted VB structures, each one be- 
ing determined by a choice of an orbital configuration (or 
orbital occupation) and a spin coupling of these orbitals. 
We consider VB structures of the form (disregarding nor- 
malization) 

active 

inactive pairs 

i'^'^) = n -lAi n («itsV - «i«it) 

active 
unpaired 

X n °jti^ac), (2) 

where dj,^ (cr =t, i) is a spin-orbital creation operator 
and |vac) is the vacuum state of second quantization. 
The VB structures arc thus made of inactive (always 
closcd-shcU) orbitals p, spin-singlet pairs of active or- 
bitals (ij), and possibly remaining unpaired spin- up ac- 
tive orbitals q. We use inactive orbitals that are either 
localized (expanded on the basis functions centered on a 
single atom), e.g. for core orbitals, or delocalized (ex- 
panded on all the basis functions of all the atoms), e.g. 
for bonds made of inactive orbitals that do not mix with 
the active orbitals. We use active orbitals that are al- 
ways localized on a single atom, and they arc typically 
identified with valence atomic (hybrid) orbitals. Note 
that Eq. ^ encompasses the case of spin-singlet pairing 
of an active orbital with itself, i.e. i = j giving sim- 
ply al^dl^ — o-l^o,]^ = ^Oi^aj^j^. The spin-coupling scheme 
based on singlet pairing used in Eq. ^ is usually referred 
to as the Heitler-London-Slater-Pauling (HLSP) scheme. 
There exist other spin-coupling schemes, but the HLSP 
scheme has the advantage of providing a clear correspon- 
dence between each VB structure and a Lewis chemical 
structure, two singlet-paired active orbitals representing 
either a bond or a lone pair. In principle, considering 
all possible pairings exhausts, for a given orbital config- 
uration, all the spin eigenstates (or spin couplings) of 
fixed quantum numbers S ~ A^unpaired/S and Ms = +S 
(A^unpaired is thc numbcr of spin- up unpaired electrons). 
In fact, considering all possible pairings leads to an over- 
complete set of spin couplings, but they can be reduced 
to a complete basis of (non-redundant) spin couplings, 
also called Rumer basis [ll,!!^,!!!]. In practice, for most 
practical applications, only a small number of chemically 
relevant VB structures are kept in the calculation. 

In practical VBSCF calculations, each VB structure 
is expanded in 2^p-''''= Slater determinants, |$/) = 
d/.plZ?^), where A^pairs is the number of pairs of dif- 
ferent active orbitals {i ^ j). The coefficients of the 
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determinants d/.^ for a given VB structure are all equal 
in absolute value. The VB structure coefficients and or- 
bital coefficients on the basis are then optimized using 
a direct generalization of the usual multi-configuration 
self-consistent-field (MCSCF) algorithm for arbitrary 
nonorthogonal orbitals [l^], or algorithms more specific 
to VB theory [H EO]. In fact, this VBSCF procedure is 
general enough to also permit optimization of GVB and 
SCVB wave functions if the active orbitals are allowed to 
delocalize. Note that GVB wave functions, in the most 
commonly used perfect-pairing approximation, are made 
of only one VB structure, while the SCVB wave functions 
include all possible spin couplings. 

As an illustration, consider the single-bonded Li2 
molecule. A typical VBSCF wave function is 

I^'vBSCf) = Ci|$cov> + C2|«'ionA) + Cal^ionB), (3) 

with a covalent VB structure obtained by spin pairing 
the 2s orbital of the ffist Li (A) with the 2s orbital of 
the second Li (B) 

l^cov) = al^AsAAssAsni 

and two ionic VB structures obtained by either spin pair- 
ing 2sA with itself, 

|$io„A) = aL.4T^l«Ai°LBt°lsB;«L^t°2s^i|vac), (5) 

or, symmetrically, by spin pairing 2s b with itself, 

l^ionB) = a\^^^a\^^^a\^^^al^^^al^^^al^^^\vac)^ (6) 

the core inactive orbitals Is a and Isb remaining always 
doubly occupied. These three VB structures correspond 
respectively to three Lewis chemical structures, Li — Li, 
~Li Li"*", and "''Li Li~, respectively. Obviously, in this 
case, spatial symmetry implies that C2 = C3 and the two 
ionic VB structures can thus be combined in a single 
"symmetry-adapted structure" , which we refer to as a 
configuration state function (CSF). In this simple case, 
the VBSCF wave function is essentially a complete ac- 
tive space (CAS) wave function with the two valence 
electrons distributed in the valence bonding ag and anti- 
bonding (T„ molecular orbitals, i.e. CAS(2,2), and would 
be strictly identical to it in a minimal basis set. In less 
trivial cases, the selected VB structures correspond only 
to a subspace of a CAS. However, if structure selection 
is correctly made, the VBSCF wave function can lead to 
observable predictions in close agreement to a CAS wave 
function. Indeed, in addition to the most chemically sig- 
nificant structures, the CAS wave function includes a lot 
of improbable covalent and (multi-)ionic structures that 
are of little importance. Thus, a VB wave function can 
represent an intelligent way of truncating a CAS expan- 
sion based on chemical grounds. 

We adopt the usual definitions of electron correlations 
used in quantum chemistry. Static correlation corre- 
sponds to the correlation included in a FVCAS wave 



function as compared with a Hartree-Fock reference, and 
dynamic correlation corresponds to the remaining corre- 
lation resulting from all the outer-valence excitations in 
a full configuration interaction expansion. The VBSCF 
method handles static correlation only. Dynamic corre- 
lation can be added in the usual way by performing con- 
figuration interaction on top of VBSCF [4l|, \^ , referred 
to as VBCL or through a perturbative treatment [4^.[43|. 
The BOVB method is an alternative way of incorporat- 
ing some dynamic correlation, keeping the wave function 
as compact as the VBSCF one. This method consists in 
allowing a different set of (inactive and active) orbitals 
in each VB structure. The orbitals can thus adapt to 
the different charge distributions in the VB structures. 
For example, for the Li2 wave function of Eqs. ([SI)-®, 
in BOVB the 2s orbitals can become more diffuse in the 
ionic structures than in the covalent structure. This extra 
degree of freedom is seen as introducing dynamic corre- 
lation, since it roughly corresponds to adding excitations 
to atomic orbitals of higher principal quantum numbers. 
Just as in the VBSCF method, inactive valence orbitals 
can be chosen to be either localized or delocalized in the 
BOVB method. Finally, in the ionic structures, a doubly- 
occupied active orbital (the 2s orbital for Li2) is some- 
times split into two independent orbitals that are spin- 
singlet coupled, which introduces additional dynamic cor- 
relation. The most flexible BOVB method is thus ob- 
tained by both splitting of active orbitals in the ionic 
structures and delocalization of inactive orbitals, which 
is referred to as split-delocalized-BOVB (SD-BOVB) in 
the VB literature. 



B. Jastrow- Valence-Bond wave functions 

We discuss now the parametrization of the two types 
of Jastrow- Valence-Bond wave functions, J- VBSCF and 
J-BOVB, that we use in QMC. The J- VBSCF wave func- 
tion is parameterized as (see Ref. 

I'f .1-vbscf) = J e^^ci\<Pi), (7) 

where J is a Jastrow operator, e'^ is an orbital "rota- 
tion" operator, and |<I>/) are CSFs, i.e. VB structures 
or linear combination of a few VB structures to accom- 
modate spatial symmetry. Each CSF is a short linear 
combination of products of spin-up and spin-down Slater 
determinants, composed generally of nonorthogonal or- 
bitals expanded on Slater basis functions. The orbital 
operator is k = ^^.j Kki Eki where Kki are the orbital 
"rotation" parameters and Eki is the singlet excitation 
operator from orbital I to orbital fc, Eki = <^If ^it + ^I^^'i' 
written with dual biorthogonal orbital creation and an- 
nihilation operators a[.^ and bia (see, e.g., Ref. 113). The 
parameters to optimize are the parameters in the Jastrow 
factor, the CSF coefficients, the orbital rotation param- 
eters and the exponents of the basis functions. 
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The J-BOVB wave function is similar except that each 
CSF is allowed to have a different set of orbitals, which 
can be written as 



:C« 



:ne 



:o~o: 



1^ 



J-BOVB 



ci e^^|$,), 



(8) 



where kj is the orbital operator for the set of orbitals 
in the CSF 1$/). In practice, Eq. ([8]) can be recast in 
the form of Eq. (O by defining the orbital operator k to 
operate on all sets of orbitals: k = J^i '^i- 

We use a Jastrow factor consisting of the exponen- 
tial of the sum of electron-nucleus, electron-electron, 
and electron-electron- nucleus terms, written as system- 
atic polynomial and Fade expansions (see also Refs. 



C. Computational details 

We start by generating a standard restricted Hartree- 
Fock (RHF) wave function with the program Gaus- 
sian [ii] and a VBSCF or BOVB wave function with 
the program XMVB [4§| . We use the core- valence triple- 
zeta quality Slater basis (CVBl) of Ema et al. [s^. For 
the XMVB calculations, each Slater function is expanded 
into 15 (core basis functions), 9 (s and p valence ba- 
sis functions), or 6 {d polarization functions) Gaussian 
functions. We use here 5 d spherical functions while 6 d 
Cartesian functions were used in Ref. 0. The standard 
RHF, VBSCF or BOVB wave function is then multiplied 
by our Jastrow factor, and QMC calculations are per- 
formed with the program CHAMP 5l| using the true 
Slater basis set rather than its Gaussian expansion. 

The Jastrow, CSF, orbital and basis exponent param- 
eters arc simultaneously optimized by minimizing the en- 
ergy plus a small percentage of the energy variance (1%) 

^"^Un VMC, us- 
~5l, m. Once 



with the linear optimization method 
ing an accelerated Metropolis algorithm 
the trial wave function has been optimized, we perform 
a DMC calculation within the short-time and fixed-node 
(FN) approximations (see, e.g., Refs. Is^ - fssl) . We use an 
imaginary time step of r = 0.01 hartree"^ in an efficient 
DMC algorithm with very small time-step errors [s^ . 

For orbital optimization, the orbitals are partitioned 
into three sets: inactive (doubly occupied in all determi- 
nants), active (occupied in some determinants and un- 
occupied in others), and virtual (unoccupied in all de- 
terminants). All inactive (including Is core) and active 
orbitals are optimized. The non-redundant excitations 
to consider are inactive — > active, inactive — ^ virtual, ac- 
tive — > virtual and active — !> active. If the action of the 
excitation E^i on the wave function is not zero but the 
reverse excitation Eik is zero, then the "orthogonality" 
condition Kik = —nu is imposed. For some active-active 
excitations, both direct and reverse excitations {E^i and 
Elk) may be allowed, and thus it makes sense for local- 
ized orbitals to remove the "orthogonality" constraint by 
treating n^i and Kjfc as independent parameters. This 



FIG. 1. Dominant bonding patterns for the C2, N2, O2, and 
F2 molecules, defining our single bonding-pattern VB wave 
functions. 



results in only a very few (if any at all) additional or- 
bital parameters for the wave functions considered here. 
We note that, when considering active-active excitations, 
redundancies between two orbital wave-function deriva- 
tives or between an orbital wave-function derivative and 
a CSF frequently occur, and must be detected and elim- 
inated. Localized orbitals do not have the point group 
symmetry of the system, so the number of orbital exci- 
tations cannot be reduced based on the non mixing of 
irreducible representations, as usually done. However, 
the number of orbital excitations is greatly reduced by 
forbidding mixing between orbitals of different localiza- 
tion classes, i.e. expanded on different subsets of basis 
functions. For BOVB wave functions, one has a differ- 
ent set of (occupied and unoccupied) orbitals for each 
CSF [g^I and only orbitals that belong to the same set 
are allowed to mix. 



III. RESULTS 



The Jastrow- Valence-Bond wave functions are tested 
on the ground states of the molecules C2 (^S+), N2 
(iS+), O2 (^S"), and F2 (^S+) at their experimental 
bond lengths of 1.2425 A, 1.0977 A, 1.2075 A, and 1.4119 
A, respectively. The core inactive Is orbitals are always 
taken as localized. For O2 and F2, we use dclocalized in- 
active TT orbitals both at the VBSCF and BOVB levels. 
This reduces the Pauli repulsion between the otherwise 
localized lone pairs. The precise composition of all wave 
functions is explained below. In addition to total ener- 
gies, we present molecular well depths all computed with 
respect to the separated atoms calculated at the FVCAS 
level. For the N2, O2, and F2 molecules, both VBSCF 
and BOVB wave functions in the limit of a complete basis 
of structures (Rumcr basis) dissociate into the FVCAS 
atomic limit, which simply corresponds for these atoms 
to single-determinant wave functions. For C2, VBSCF in 
the limit of a complete basis of structures also dissociates 
into the FVCAS atomic limit, which corresponds for the 
carbon atom in its triplet ground state to a CAS(2,2) 
wave function, since the carbon lone pair can doubly oc- 
cupy either a 2s or a 2p orbital. With this definition, 
the well depths actually allow a comparison of the abso- 
lute quality of the different molecular wave functions at 
equilibrium since the atomic reference is the same. 
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TABLE I. Relative contributions of the different types of elec- 
tron correlation to the equilibrium well depths for C2, N2, O2, 
and F2. Starting from the RHF well depth, the "left-right" 
static correlation is defined as the contribution recovered by a 
single-bonding pattern VBSCF calculation, the other parts of 
static correlation as the additional contributions recovered by 
a FVCAS MCSCF calculation, and the dynamic correlation 
as the remaining contribution necessary to reach the exact 
well depth. 

Static correlation Dynamic correlation 



left-right otlicr 



C2 


77% 


18% 


5% 


N2 


70% 


14% 


16% 


02 


60% 


8% 


32% 


F2 


60% 


8% 


32% 



A. Single bonding-pattern wave functions 

We first consider wave functions describing what we 
call a single bonding pattern. A bonding pattern is a 
compact Lewis-style picture of the molecular electronic 
structure, but which is made of several (covalent and 
ionic) VB structures. The dominant bonding patterns of 
C2, N2, O2, and F2 arc depicted in Fig. [T] The conven- 
tions of the drawings arc the following. Each single line 
connecting two atoms represents a cr or tt two-electron 
bond. In our wave functions, each two-electron bond is 
generally composed of one covalent and two ionic com- 
ponents. Hence, F2 (single bond) is described by a wave 
function containing 3 VB structures, whereas N2 (triple 
bond) leads to a wave function with 3^ = 27 VB struc- 
tures. Each pair of dots on a single atom represents a 
lone pair, which is described by two opposite-spin elec- 
trons occupying the same orbital. 

The single dots connected by a curve which appears 
for C2 represent a (formally covalent) singlet coupling 
between two electrons occupying two different orbitals, 
one located on each carbon in that case. The two corre- 
sponding coupled orbitals are sp hybrids pointing in op- 
posite direction, so that they overlap and thus interact 
very weakly. Hence, this corresponds to a singlet birad- 
ical state, the coupling between the two single electrons 
being only formal to ensure the proper spi n symmetry. 
Very recent VB calculations by Su et al. |62| have shown 
this biradical triple-bond situation to be the dominant 
bonding pattern for C2 at equilibrium distance. In or- 
der to keep the number of VB structures low, only the 
single-ionic structures, as well as the double-ionic struc- 
tures with opposite charges on each carbon, are included 
in our wave function. This leads to a wave function with 
21 VB structures. 

The O2 molecule is considered in its spin-triplet (para- 
magnetic) ground state. The picture in Fig. [1] corre- 
sponds to the textbook molecular-orbital diagram for this 
molecule, displaying a u two-electron bond (single line 
connecting the two oxygens), and two tt three-electron 
half-bonds, each one represented by a three-dot symbol. 
In our wave functions, the a bond is composed of one co- 



valent and two ionic components, and each three-electron 
bond is composed of the two resonating structures typ- 
ical of this kind of bond [3l|, [s^] ■ This leads to a wave 
function with 12 VB structures. 

In each wave function, structures equivalent by sym- 
metry are further grouped together into CSFs. For each 
dimer, the VBSCF and BOVB wave functions have ex- 
actly the same number of CSFs, the only difference being 
that different sets of orbitals are used for the different 
CSFs of the BOVB wave functions. 

A rough analysis of the relative contributions of the 
different types of electron correlation to the equilibrium 
well depths for the four molecules is given in TableU) This 
analysis leads us to consider separately the case of the O2 
and F2 molecules which have strong dynamic correlation, 
and the case of the C2 and N2 molecules which have 
strong static correlation. 



1. Systems with strong dynamic correlation: O2 and F2 

The total energies and well depths of O2 and F2 com- 
puted with standard VB methods and QMC methods are 
displayed in Table [Hi We start by analyzing the standard 
VB results. The difference between the RHF and VBSCF 
energies can be taken as the definition of the "left-right" 
part of static correlation. In VB theory, the left-right 
correlation is viewed as the energy gained by restoration 
of the correct balance between covalent and ionic con- 
tributions to the bond, starting from the incorrect 50% 
covalent - 50% ionic RHF description. The contributions 
of left-right correlation to the well depths appear to be 
very large here, 59.3 and 45.1 kcal/mol for O2 and F2, 
respectively. A significant part of dynamic correlation is 
retrieved by going from VBSCF to BOVB, improving well 
depths by 33.2 and 18.9 kcal/mol for O2 and F2, respec- 
tively. Another 3.9 kcal/mol improvement is obtained for 
F2 by splitting the doubly-occupied active orbital in the 
ionic structures, which brings radial correlation to the 
active electron pair. The agreement between the BOVB 
and exact well depths is quite reasonable considering the 
compactness of these wave functions; only 12 VB struc- 
tures (combined into 5 CSFs) for O2 and 3 VB structures 
(combined into 2 CSFs) for F2. 

We now discuss the VMC results. The Jastrow factor 
keeps electrons apart and describes part of dynamic and 
static correlation. It is a particularly efficient way of in- 
cluding dynamic correlation, that otherwise requires in 
configuration interaction calculations a very large num- 
ber of highly excited determinants. It also includes some 
static left-right correlation since it greatly reduces the 
probability of ionic configurations |6|- As expected, a 
sharp decrease of the total energies is observed in VMC, 
indicating that a large amount of the total correlation 
energy is recovered by the Jastrow factor. However, for 
well depths, and most chemical properties, we are in- 
stead interested in a balanced description of the differen- 
tial correlation energy in the molecule and the separated 
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TABLE II. Total energies (in hartree) and well depths (in kcal/mol) for the ground states of O2 and F2 at their experimental 
bond lengths calculated in RHF, VBSCF and BOVB, and in VMC and DMC using corresponding Jastrow- Valence-Bond 
wave functions with all the parameters (Jastrow, CSF coefficients, orbitals and basis exponents) fully optimized in VMC. For 
comparison, DMC results for FVCAS wave functions from a previous study are also shown. 

Total energies Well depths Total energies Well depths 



RHF 


-149 


6567 


22.4 


-198 


7608 


-36.2 


VBSCF 


-149 


7512 


81.7 


-198 


8327 


8.9 


BOVB'= 


-149 


8042'' 


114.9'* 


-198 


8627 


27.8 


BOVB-splif^ 








-198 


8689 


31.7 


VMC J-RHF 


-150 


2287(5) 


99.0(4) 


-199 


4204(5) 


12.4(4 


VMC J- VBSCF 


-150 


2414(5) 


107.0(4) 


-199 


4400(5) 


24.7(4 


VMC J-BOVB'' 


-150 


2496(5) 


112.1(4) 


-199 


4456(5) 


28.3(4 


VMC J-BOVB-split = 








-199 


4494(5) 


30.6(4 


DMC J-RHF 


-150 


2861(2) 


114.7(2) 


-199 


4833(2) 


29.8(2 


DMC J-VBSCF 


-150 


2889(2) 


116.5(2) 


-199 


4900(2) 


34.0(2 


DMC J-BOVB'= 


-150 


2922(2) 


118.5(2) 


-199 


4922(2) 


35.4(2 


DMC J-BOVB-split= 








-199 


4962(2) 


37.9(2 


DMC J-FVCAS" 


-150 


2944(1) 


119.6(1) 


-199 


4970(1) 


37.9(1 


Estimated cxaet^ 


-150 


3274 


120.9(<1) 


-199 


5304 


39.0(1 



" Ref.JS 
Rcf. [gj. 

BOVB with dolocalized vr orbitals, which is referred to as tt-D-BOVB in the VB literature. 

Starting from an optimized BOVB wave function with all localized orbitals, delocalizcd tt orbitals are allowed and optimized but 
keeping all the other orbitals frozen. Optimization of all the orbitals with delocalizcd tt orbitals was not stable. This problem was not 
encountered in the VMC optimization. 

"BOVB-split" (usually referred to as SD-BOVB in the VB literature): in each of the ionic structures, the doubly occupied active 
orbital is split into two singlet-coupled singly occupied orbitals. 



atoms. The well depths improve considerably going from 
RHF to J-RHF, by 76.6 and 48.6 kcal/mol for O2 and F2, 
respectively, confirming that the Jastrow factor recovers 
as well a significant part of left-right static correlation. 
Consequently, going from J-RHF to J-VBSCF wave func- 
tions improves the VMC well depths by only 8.0 and 11.7 
kcal/mol for O2 and F2, respectively, much less that is 
gained by going from standard RHF to standard VBSCF. 
The J-VBSCF VMC well depths remain in relatively poor 
agreement with the exact values, underestimated by 13.9 
and 14.3 kcal/mol for O2 and F2, respectively. It is ex- 
pected that these errors are mainly due to some missing 
dynamic correlation, rather than static correlation corre- 
sponding to some missing VB structures. This point will 
be checked in Section III.B. 

Additional correlation can be introduced in VMC by 
replacing the VBSCF determinant expansion by the more 
flexible BOVB form. By using a different set of orbitals 
in each VB structure, the BOVB wave function allows 
the orbitals to be more diffuse on the negatively charged 
atom, and more compact on the positively charged atom 
of the ionic structures, compared to the neutral covalent 
structure. This "breathing-orbital effect" can be seen 
as a manifestation of dynamic correlation. It is clear 
that such breathing-orbital effects are partly redundant 
with the correlations introduced by the Jastrow factor. 
In practice, it is observed that with the limited form of 
the Jastrow factor employed, the J-BOVB wave func- 
tions improve a bit upon J-VBSCF wave functions: the 



VMC well depths are bettered by 5.1 and 2.5 kcal/mol 
for O2 and F2, respectively. As expected, these improve- 
ments are much smaller than those obtained when go- 
ing from standard VBSCF to BOVB, showing that a 
large part of dynamic correlation is already described by 
the Jastrow factor. The agreement between the VMC 
J-BOVB well depths and the exact values is still not 
good, with an underestimation of 8.8 kcal/mol for O2, 
and 10.7 kcal/mol (8.4 kcal/mol with BOVB-split) for 
F2, respectively. Even more disappointing, the VMC J- 
BOVB well depth is actually worse than the standard 
BOVB one for O2, and the same is true for BOVB-split 
for F2. This indicates that our Jastrow factor describes 
correlation effects in the atom more accurately than in 
the (less symmetric) diatomic molecule. By incorporat- 
ing more correlation in the separated atoms than in the 
molecule, the Jastrow factor thus makes the BOVB well 
depths less accurate. In fact, it is significant that for 
F2 splitting the ionic electron pairs improves the VMC 
well depth by 2.3 kcal/mol, indicating that this radial 
correlation effect in the molecule was not well described 
by the Jastrow factor. The more flexible multi- Jastrow 
approach of Ref. [l^ could solve this problem by better 
handling the delicate balance between atomic and molec- 
ular correlation effects [6^ . 

Let us now consider the fixed-node DMC results. The 
total energies and the well depths reflect the quality of 
the nodes of the different wave functions. Both total 
energies and well depths are much improved compared 
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to the VMC results. J-RHF wave functions give DMC 
well depths of moderate accuracy, underestimated by 6.2 
and 9.2 kcal/mol for O2 and F2, respectively. J-VBSCF 
wave functions improves DMC well depths but by rather 
limited amounts, 1.8 and 4.2 kcal/mol for O2 and F2, 
respectively, showing that most of left-right correlation 
was already recovered with the J-RHF nodes. Using J- 
BOVB wave functions allow further gains in accuracy, 
leading to results close to the ones obtained with much 
less compact J-FVCAS wave functions. The best DMC 
J-BOVB calculations underestimate the well depths by 
only 2.4 and 1.1 kcal/mol for O2 and F2, respectively. 



2. Systems with strong static correlation: C2 and N2 

Let us now consider C2 and N2, two molecules hav- 
ing strong static correlation even at equilibrium dis- 
tance. RHF is expected to be particularly poor for these 
molecules, because of the triple bond nature of these sys- 
tems. In a RHF wave function, each two-electron bond 
is constrained to be 50% covalent - 50% ionic, and in ad- 
dition all (multi-)ionic forms are constrained to have the 
same weights, which is unphysical. In particular, bi-ionic 
structures with two positive charges on one atom and two 
negative charges on the other will have the same weight as 
neutral bi-ionic structures (one positive and one negative 
charge on each atom) . Similarly, tri- ionic structures with 
three (positive or negative) charges of the same sign on 
each atom will be forced to have the same weight as the 
mono-charged-atom tri-ionic structures. Note that GVB 
wave functions also have these constraints on the weights 
of multi-ionic structures, this form of VB wave function 
being only able to take into account intrapair left-right 
correlation (covalent vs. ionic balance within only one 
two-electron bond) but not interpair correlation. On the 
other hand, a VBSCF wave function can include both in- 
trapair and interpair left-right correlation. More impor- 
tantly, for C2, the singlet biradical character is expected 
to generate additional static correlation beyond left-right 
correlation, which will require extra bonding patterns. 

The total energies and well depths of these molecules 
are displayed in Table IIIII As expected, RHF strongly 
underestimates the well depths (note that the negative 
well depth obtained for C2 is in fact due to the choice of 
the FVCAS reference for the separate atoms). VBSCF 
appears to retrieve a large part of the missing correlation, 
as it improves the well depths by 120.6 and 78.8 kcal/mol 
for C2 and N2, respectively, as compared with RHF val- 
ues. The VBSCF well depths are however still underes- 
timated by more than 30 kcal/mol for both molecules. It 
was impossible to carry out standard BOVB calculations 
to assess dynamic correlation effects at the VB level, as 
the number of VB structures (21 for C2 and 27 for N2) 
is too large to allow such calculations. Last, we checked 
that GVB (not shown in Table IIIip greatly underesti- 
mates the well depths as expected, giving 84.0 and 164.9 
kcal/mol for C2 and N2, respectively, values which are 



about 30 kcal/mol worse than those from VBSCF. 

As for O2 and F2, addition of a Jastrow factor in 
VMC considerably improves the results. This improve- 
ment is dramatic for single-determinant wave functions, 
the J-RHF weU depths being bettered by 110.7 and 88.4 
kcal/mol for C2 and N2, respectively, as compared with 
the RHF ones. The improvement is much less impor- 
tant but still significant for the VBSCF wave functions, 
with J-VBSCF well depths improved by 21.6 and 22.6 
kcal/mol for C2 and N2, respectively, compared to the 
VBSCF ones. As for O2 and F2, these findings show 
that the Jastrow factor can retrieve a significant amount 
of left-right correlation. The accuracy of the VMC J- 
VBSCF calculations is similar to what has been found 
for O2 and F2, with well depths underestimated by 15.1 
and 11.3 kcal/mol for C2 and N2, respectively. 

The DMC resuhs follow the trends found in VMC, but 
the differences between C2 and N2 show up more clearly. 
The N2 well depth is underestimated by 7.3 kcal/mol 
with a J-RHF wave function, which is actually in the 
range of what was observed for O2 and F2. A J-VBSCF 
wave function leads to a DMC well depth close to chemi- 
cal accuracy, indicating that the quality of the J-VBSCF 
nodes are good enough get good energy differences. The 
situation is quite different for C2. The J-RHF DMC well 
depth is underestimated by 26.7 kcal/mol, far from the 
accuracy reached for the other molecules. Using a J- 
VBSCF wave function greatly improves the DMC well 
depth, but with an error of 8.6 kcal/mol, it is still far 
from chemical accuracy. For O2 and F2; a further im- 
provement is possible by using J-BOVB wave functions, 
but for C2 even if it were possible to use a J-BOVB wave 
function a significant improvement is not expected. Be- 
cause of the multiradical character of C2 , we expect the 
missing correlation to still be of static nature, rather than 
dynamic. The route for improvement is to go beyond a 
single bonding-pattern description. 



B. Multiple bonding-pattern wave functions 

A straightforward way of improving our J-VBSCF 
wave functions, particularly for the C2 molecule, is to 
increase the number of VB structures included in the 
wave function. Adding VB structures should be done in 
a systematic but selective way, as the total number of 
non-redundant VB structures (for all orbital configura- 
tions and spin couplings) can be very large, e.g. up to 
1764 structures for C2 (for the formula giving the num- 
ber of VB structures in a complete (Rumer) basis, see 
Ref. [U). A route to extend the VB description beyond 
a single bonding pattern is to systematically derive pos- 
sible alternative bonding patterns, and classify them in 
terms of stability following Lewis-type rules. In particu- 
lar, the most bonded patterns will be preferred over the 
less bonded ones, a a bond will be preferred over a tt 
one, and a three-electron (or one-electron) bond will be 
counted as half a bond. Following these rules of thumb. 
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TABLE III. Total energies (in liartree) and well depths (in kcal/mol) for the ground states of C2 and N2 at their experimental 
bond lengths calculated in RHF and VBSCF, and in VMC and DMC using corresponding Jastrow- Valence-Bond wave functions 
with all the parameters (Jastrow, CSF coefficients, orbitals and basis exponents) fully optimized in VMC. For comparison, DMC 
results for FVCAS wave functions from a previous study are also shown. 



C2 N2 

Total energies Well depths Total energies Well depths 



RHF 


-75.4014 


-8.8 


-108.9856 


115.8 


VBSCF 


-75.5936 


111.8 


-109.1112 


194.6 


VMC J-RHF 


-75.8132(5) 


101.9(4) 


-109.4516(5) 


204.2(4) 


VMC J-VBSCF 


-75.8635(5) 


133.4(4) 


-109.4723(5) 


217.2(4) 


DMC J-RHF 


-75.8667(2) 


121.8(2) 


-109.5039(2) 


221.2(2) 


DMC J-VBSCF 


-75.8954(2) 


139.9(2) 


-109.5119(2) 


226.3(2) 


DMC J-FVCAS" 


-75.9106(1) 


149.5(1) 


-109.5206(1) 


231.5(1) 


Estimated exact** 


-75.9265 


148.5(5) 


-109.5427 


228.5(<1 



" Ref.ll 
>> Rcf . [gl 



a description extended to multiple bonding patterns is 
proposed in Fig. [2]for the C2, N2, O2, and F2 molecules. 

The main bonding pattern for C2 is of acetylenic type, 
displaying one a and two tt bonds, as well as two singlet- 
coupled (biradical) electrons. It is natural to postu- 
late for this molecule a secondary bonding pattern of 
ethylenic type, with one a and one tt bond, thus leaving 
one pair of singlet-coupled electrons on each atom. An- 
other doubly bonded pattern could be postulated for C2, 
with two TT bonds. However, by applying our simple rules, 
it is expected to be less stable than the ethylenic pattern, 
as a bonds are usually more stable than tt ones. In a re- 
cent VB study, this 7r-doubly bonded pattern has indeed 
been found to be of very little importance for the ground 
state at equilibrium distance (G^I , so it will not be consid- 
ered here. Considering a C2 description with two bonding 
patterns (acetylenic and ethylenic) , and including the two 
covalent and all mono-ionic and neutral bi-ionic struc- 
tures, leads to a set of 35 VB structures, which can be 
combined into only 12 CSFs when all the structures that 
must have equal weights by symmetry are combined. All 
calculated well depths arc gathered in Table IIVI The 
first three lines just repeat, for ease of comparison, the 
well depths obtained using single bonding-pattern (J- 
)VBSCF wave functions, previously shown in Tables HIl 
and mil The next three lines display the well depths ob- 
tained using multiple bonding-patterns (J-)VBSCF wave 
functions, and the last three lines some reference data. At 
the standard VBSCF level already, the well depth of C2 
is significantly improved, by 15.4 kcal/mol, when going 
from the single bonding-pattern to the multiple bonding- 
pattern wave function. Interestingly, a large part of this 
gain is retained at the VMC J-VBSCF level, with a 9 
kcal/mol improvement. This suggests that the Jastrow 
function is not able to significantly recover the part of 
static correlation that is included when going from sin- 
gle bonding pattern to multiple bonding patterns. This 
observation partly transfers to the DMC results as well, 
with a significant 4.9 kcal/mol improvement of the well 



depth. In conclusion, the multiple bonding-pattern J- 
VBSCF wave function for C2 gives both VMC and DMC 
well depths in reasonable agreement with estimated exact 
values (errors of 6.1 and 3.7 kcal/mol, respectively), while 
still being a very compact wave function (only 12 CSFs 
compared to 165 CSFs for the FVCAS wave function). 
The remaining error could be progressively eliminated by 
adding VB structures corresponding to additional bond- 
ing patterns, starting with the pattern with two tt bonds 
and no a bond, and then including patterns with only 
one (cr or tt) bond. It is expected that a J-VBSCF wave 
function including all the main VB structures belonging 
to these bonding patterns would give almost exact disso- 
ciation energies, as has recently been shown to be the case 
at the VBCI level [s^, while still having a significantly 
more compact wave function than the FVCAS one. 

Following the same reasoning as for C2, a secondary 
bonding pattern displaying one a and one tt bond can be 
postulated for the N2 molecule, corresponding to picture 
b in Fig. [51 This bonding pattern could also be obtained 
by starting from the common Lewis structure of spin- 
singlet O2, and doubly ionizing one of the two lone pairs. 
For N2 , a TT-doubly bonded pattern is likely to be of even 
less importance than for C2, as this would result in an 
N atom with a completely filled a space, i.e. two (highly 
repulsive) a lone pairs facing each other, so we rule it 
out. Last, singly bonded patterns with either a cr or a tt 
bond, corresponding to picture c, are also included, by 
analogy with the secondary bonding pattern of C2. For 
all these extra patterns, the covalent and all the neu- 
tral and mono-chargcd-atom ionic VB structures arc in- 
cluded. This gives 27 extra VB structures, and so a total 
of 54 structures (with the main bonding pattern) which 
can be combined into 24 CSFs. This multiple bonding- 
pattern wave function improves the VBSCF well depth of 
N2 by 15.4 kcal/mol in comparison to the single bonding- 
pattern calculation, while still being a substantially more 
compact wave function than the FVCAS one (24 vs. 107 
CSFs). The improvements in VMC and DMC are less 
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FIG. 2. Most important resonant bonding patterns for the 
C2, N2, O2, and F2 molecules, defining our multiple bonding- 
pattern VB wave functions. 



TABLE IV. Well depths (in kcal/mol) for the ground states of 
C2, N2, O2 and F2 at their experimental bond lengths calcu- 
lated in VBSCF, VMC J-VBSCF and DMC J-VBSCF using 
single and multiple bonding-pattern wave functions. All the 
parameters of the J-VBSCF wave functions (Jastrow, CSF 
coefficients, orbitals and basis exponents) have been fully op- 
timized in VMC. For comparison, results for FVCAS wave 
functions from a previous study are also shown. 



C2 



N2 



O2 



Single bonding pattern 
111.8 194.6 81.7 8.9 
133.4(4) 217.2(4) 107.0(4) 24.7(4) 
139.9(2) 226.3(2) 116.5(2) 34.0(2) 



VBSCF 

VMC J-VBSCF 
DMC J-VBSCF 



VBSCF 

VMC J-VBSCF 
DMC J-VBSCF 



MCSCF FVCAS" 
VMC J-FVCAS° 
DMC J-FVCAS" 
Estimated cxact^ 



Multiple bonding patterns 
127.2 203.4 84.1 9.6 
142.4(4) 220.2(4) 108.2(4) 26.9(4) 
144.8(2) 228.1(2) 117.5(2) 35.7(2) 

References 

140.8 210.0 89.9 15.3 

146.9(1) 225.5(2) 108.7(2) 27.4(2) 

149.5(1) 231.5(1) 119.6(1) 37.9(1) 

148.5(5) 228.5(<1) 120.9(<1) 39.0(1) 



" Ref.la 
Ref . m 



pronounced, though not negligible, with a 2-3 kcal/mol 
gain. The DMC J-VBSCF well depth is very close to the 
estimated exact value, indicating the good quality of the 
nodes of the multiple bonding-pattern wave function. 

For O2, the main bonding pattern corresponds to a 
double-bond situation, with a a bond and two tt half 
bonds. There are not many possibilities for deriving al- 
ternative patterns that are compatible with the spatial 
and spin-triplet symmetry of the ground state, and in 
particular it is impossible to preserve the a bond. Pat- 



tern b keeps the two three-electron bonds, while pattern c 
has one three-electron and one one-electron tt half-bonds 
and a completely filled a space on each atom. There 
are a total number of 12 VB structures emerging from 
patterns b and c, and when added to the previous set 
we arrive at a total number of 24 structures (10 CSFs). 
Obviously, even if these extra patterns can be associated 
with a formal bond order of one, they are not expected 
to be of much importance in the ground state. This ex- 
pectation is already confirmed at the standard VBSCF 
level, as the well depth is improved by only 2.4 kcal/mol, 
much less than what was observed for C2 and N2. Quite 
logically, the gains are even smaller in VMC and DMC, 
with about 1 kcal/mol improvement only. It thus seems 
that for O2 the single bonding-pattern wave function al- 
ready contains most of the static correlation, and that 
further gain in accuracy requires instead improvement 
of the description of dynamic correlation, as previously 
found when going from J-VBSCF to J-BOVB wave func- 
tions. 

For F2, alternatives to the usual cr bond situation are 
limited too. A tt bond situation can be postulated, how- 
ever it comes with a (highly repulsive) completely filled a 
space on each atom. A situation with a a three-electron 
bond can be envisaged (picture b) as it corresponds to the 
(quite stable) ground state of the difluorine anion, but it 
comes here with a singly-ionized a lone pair. Consider- 
ing all these bonding patterns generates a set of 13 VB 
structures (5 CSFs). As expected, the multiple bonding- 
pattern wave function gives very small improvements. As 
for O2, a much more significant improvement was previ- 
ously obtained for F2 by using a BOVB wave function. 



C. Effect of orbital optimization 

Although for some systems reoptimizing the orbitals in 
QMC in the presence of the Jastrow factor leads to only 
small improvements, it is known that it can be very im- 
portant for other systems, e.g. for calculating certain ex- 
citation energies |65l - [67| . To examine the effect of orbital 
optimization, we show in Table |V] the VMC and DMC 
well depths of the four molecules for J-RHF and single 
bonding-pattern J-VBSCF wave functions with or with- 
out VMC optimization of the orbital coefficients and the 
basis exponents. In each case, of course, the same level of 
optimization has been used for the atoms. At the VMC 
level, optimization of the orbital coefficients and the ba- 
sis exponents improves the well depths by the order of 
10 kcal/mol for J-RHF and J-VBSCF wave functions. In 
DMC, the gain from optimization of the orbital coeffi- 
cients and the basis exponents is smaller but still signif- 
icant, ranging from about 2 to 7 kcal/mol. This shows 
that orbital optimization in QMC has a significant im- 
pact of the nodes of the wave functions. These improve- 
ments in VMC and DMC are primarily due to the opti- 
mization of the orbital coefficients, the gains from reopti- 
mization of the basis exponents being generally smaller. 
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TABLE V. Well depths (in kcal/mol) for the ground states of C2, N2, O2 and F2 at their experimental bond lengths calculated 
in VMC and DMC using single bonding-pattern J-VBSCF wave functions with partial or complete VMC optimization of the 
parameters: Jastrow (J), CSF coefficients (c), orbital coefficients (o) and basis exponents (e). 





optimized param. 


C2 


N2 


O2 


F2 


VMC J-RHF 


J 


86.4(4) 


191.0(4) 


84.5(4) 


-0.2(4) 


VMC J-RHF 


J+o 


100.9(4) 


202.7(4) 


91.0(4) 


11.7(4) 


VMC J-RHF 


J+0+0 


101.9(4) 


204.2(4) 


99.0(4) 


12.4(4) 


VMC J-VBSCF 


J+c 


125.6(4) 


207.6(4) 


95.8(4) 


13.4(4) 


VMC J-VBSCF 


J+c-l-o 


129.5(4) 


213.8(4) 


100.0(4) 


24.0(4) 


VMC J-VBSCF 


J+c+o-fe 


133.4(4) 


217.2(4) 


107.0(4) 


24.7(4) 


DMC J-RHF 


J 


114.5(2) 


218.8(2) 


110.9(2) 


25.0(2) 


DMC J-RHF 


J+o 


122.1(2) 


219.7(2) 


113.0(2) 


30.8(2) 


DMC J-RHF 


J+o+c 


121.8(2) 


221.2(2) 


114.7(2) 


29.8(2) 


DMC J-VBSCF 


J+c 


135.5(2) 


222.7(2) 


113.8(2) 


30.1(2) 


DMC J-VBSCF 


J+c+o 


138.5(2) 


224.4(2) 


113.6(2) 


33.4(2) 


DMC J-VBSCF 


J+c+o+e 


139.9(2) 


226.3(2) 


116.5(2) 


34.0(2) 


Estimated cxact^ 




148.5(5) 


228.5(<1) 


120.9(<1) 


39.0(1) 



TABLE VI. Numbers of CSFs (A^csf), numbers of unique spin-up -I- spin-down determinants (A'dct) in ground-state wave 
functions and numbers of unique spin- up -I- spin-down singly excited determinants (A'^cO for orbital optimization for the single 
bonding-pattern (SBP) and multiple bonding-pattern (MBP) VBSCF and BOVB wave functions employed in this work. Values 
for FVCAS wave functions are also shown for comparison. 







C2 






N2 






O2 






F2 






NcSF 


JVdot 




NcSF 


Ndot 


"dot 


NcSF 




-"^dot 


NcSF 


Wdct 


-"^dot 


SBP VBSCF 


a 


32 


1488 


10 


16 


1280 


5 


10 


848 


2 


4 


232 


SBP BOVB 














5 


32 


2144 


2 


20 


1428 


MBP VBSCF 


12 


64 


2992 


24 


72 


5152 


10 


28 


2282 


5 


16 


1572 


FVCAS 


165 


140 


6120 


107 


112 


5864 


30 


52 


2548 


8 


16 


984 



except in the case of 02- 



D. Computational cost 

Finally, we discuss the computational cost of our 
Jastrow- Valence-Bond wave functions. In QMC calcula- 
tions, using multideterminant wave functions, the major 
bottleneck is the evaluation of Slater determinants along 
with their first- and second-order derivatives with respect 
to electron coordinates, so that the computational cost is 
usually directly determined by the number of these deter- 
minants. In addition, orbital optimization in VMC using 
a super-configuration-interaction- type algorithm [3,01 re- 
quires evaluation of numerous singly excited determi- 
nants, which is often the bottleneck of the optimization 
calculations. It is thus desirable to have compact wave 
functions with small numbers of determinants. In prac- 
tice, each determinant is actually broken into the prod- 
uct of a spin-up and a spin-down determinant, and only 
the list of unique spin-up and spin-down determinants is 
computed. In Table IVTl we have reported the numbers of 
CSFs, numbers of unique spin-up -I- spin-down determi- 
nants, and the numbers of unique spin-up -|- spin-down 
singly excited determinants for the VBSCF and BOVB 
wave functions employed in this work. It is seen that 
there are about 4-7 times less determinants to evaluate 
for the single bonding-pattern VBSCF wave functions. 



in comparison to the FVCAS wave functions. For orbital 
optimization, these single bonding-pattern VBSCF wave 
functions also generate about 3-5 times less singly excited 
determinants. On the other hand, single bonding-pattern 
BOVB wave functions do not appear as computationally 
advantageous in comparison to FVCAS wave functions. 
For O2, the BOVB wave function has only a few less 
determinants, and for F2 the BOVB wave function in 
fact has more determinants than the FVCAS wave func- 
tion. Multiple bonding-pattern VBSCF wave functions 
give less determinants than FVCAS wave functions for 
C2, N2; and O2, but the reduction is rather modest with 
at most about a factor of 2 less determinants for C2 . For 
F2, the multiple bonding-pattern VBSCF wave function 
actually generates more singly excited determinants than 
the FVCAS wave function. This happens because for 
these small F2 wave functions, it is more advantageous 
to restrict the number of single excitations based on the 
spatial symmetry of the molecular orbitals, as done in the 
FVCAS calculations, rather than based on a localization 
criterion of the VBSCF orbitals. 



IV. CONCLUSIONS 

We have performed VMC and DMC calculations us- 
ing compact Jastrow- Valence-Bond wave functions based 
on strictly localized active orbitals. The J-VBSCF wave 
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functions are based on a chemically intuitive description 
of bonding, including either a single bonding pattern 
(Fig. [T|) or multiple bonding patterns (Fig. [2]). The J- 
BOVB wave functions introduce additional correlation 
through orbital relaxation. 

For N2, O2, and F2, the single bonding-pattern J- 
VBSCF wave functions yield DMC equilibrium well 
depths with errors ranging from about 2 to 5 kcal/mol. 
This is quite a reasonable accmracy considering the com- 
pactness of the wave functions. For C2, a system with 
very strong static correlation, a two-bonding-pattern J- 
VBSCF wave function is necessary to reach a similar ac- 
curacy. Optimization of all the wave-function parameters 
in VMC, and in particular the orbitals, is important to 
reach this accuracy. As regards computational cost, these 
VBSCF wave functions contain significantly less Slater 
determinants than FVCAS expansions, and are therefore 
more efficient to evaluate. The J-BOVB wave functions 
are much more costly to evaluate than the J- VBSCF wave 
functions, and give only moderate improvement of well 
depths. 

The VB ansatz provides a general way to design 
compact wave functions that capture only the essential 
physics of the chemical bond, by including only the most 
chemically relevant structures. These wave functions 



could be further improved, by including more and more 
structures belonging to secondary bonding patterns, de- 
pending on the accuracy needed. 

For first-row diatomic molecules this approach is of 
limited interest since it is easily possible to use FVCAS 
expansions. However, for molecules requiring larger ac- 
tive spaces, very compact valence bond wave functions 
can still be designed, and this becomes a great advan- 
tage over molecular orbital based multi-determinant ap- 
proaches. Also, for large systems with regions of inter- 
est localized on a few atoms, which is usually the case 
in chemical reactions, strictly localized orbitals enable 
both a reduction of the active space and a reduction of 
the orbital optimization cost, as compared to delocalized 
molecular orbital approaches. 
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